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een) ON OJ Toe Cc: 
ABSTRACT 


This thesis considers the influence of four secondary effects 
on the analysis of a simple ell-shaped piping structure when it is 
uniformly heated. The four effects are: (a) axial deformation, 
(b) shearing deformation, (c) beam-column effect, and (d) dif- 
ference between arc and chord. They are usually neglected in a 
piping flexibility analysis. 

In a paper of several years ago, J. E. Brock developed a theory 
for this investigation but <i dame tapyeteldt any numerical results or 
conclusions. The present thesis reviews Brock's theory and corrects 
some errors in coefficients. It then describes a digital computer 
program for IBM System 360 which performs the corresponding calculations. 
Finally, it draws the new and quite significant conclusion that con- 
ventional piping stress analysis, which neglects these secondary effects, 
can result in gross errors in evaluating stresses in piping config- 


urations which are likely to be encountered in piping practice. 
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SECTION I 


INTRODUCTION 


l-l1 Background 

In 1964 Professor John E. Brock of the Mechanical Engineering 
Department of the Naval Postgraduate School, Monterey, California, 
published a paper entitled, "Some Secondary Effects in a Simple 
Piping Structure Under feet 

The secondary effects referred to are axial deformation, shear, 
beam column action, and difference between arc and chord. The 
structural element considered is an ell-shaped piping configuration 
with a square corner. 

In the normal analysis of such a structure, the above mentioned 
secondary effects are neglected, being regarded as insignificant in 
relieving the stresses caused by thermal expansion. Only the 
principal sources of relief, bending and torsion, are usually consid- 
ered, 

Professor Brock undertook to make a complete flexibility analysis 
of this piping configuration including all of the four secondary 
effects. The 1964 paper reports his results. It also reports that 
a computer program had been written to investigate this analysis 
and the resulting system of equations. The program had been written 
but was not developed sufficiently to thoroughly investigate the 


relationships involved. 


John E. Brock, "Some Secondary Effects In A Simple Piping Structure 


Under Heating," Journal of Applied Mechanics, Vol 31, March 1964, 
pp. 88-90. 


1-2 Objectives of Work 

The object of this thesis was to write a computer program for 
a large digital computer that would be capable of performing the 
calculations required by Professor Brock's theoretical development 
and one that would allow application of the theory to any engineering 
problem that involves a piping structure of the given configuration. 
Further, it attempts to determine under what conditions these 
secondary effects may become significant in calculation. 
1-3 Assumptions and Limitations 

The investigations that are reported herein are based on an ell- 
shaped piping configuration with a square corner. No other shape is 
considered. Also, the configuration is assumed to have deformation 
only in the plane of the legs a and b. 

It is further assumed that deformations and displacements are 
small compared to the longitudinal dimensions of the pipe, except 
that a first approximation is included for the difference between 


arc length and chord length of the deformed pipes. 
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1-4 


Notation 


TABLE I 
NOMENCLATURE 
length of straight pipe elements 
cross sectional area of pipe material 
cross sectional area of pipe contents 
unit thermal strain 
Young's modulus of elasticity 
structural forces in ell-configuration 
coefficients defined in equations (16)-(19) 
shearing modulus of elasticity 
(1-2u) p%/e (i- p*) = pressure coefficient 
moment of inertia of pipe cross section 
length of pipe 
bending moments 
internal pressure 
lateral force 


axial compressive force 


y 


Ue 
(-us) 

QL? /EI 

axial coordinate 

horizontal deflection of joint 
lateral coordinate 

vertical deflection of joint 


coefficients (see Table II) 


al 


TABLE I (cont.) 


C, Re coefficients defined in equations (15A) 
and (15B) 
¥ CP/AG = effective shearing strain 


tip axial extension 
tip lateral deflection 
ly = conveniently small positive numbers 
shear distribution factor (see equation (9)) 
2 
G EL/AGL 


rotation of joint 


y On wr 


difference between arc and chord (first 
approximation) 


MAM, Mz coefficients in equations (9), (10), and (11) 


pt 


Poisson's ratio 
(inside diameter of pipe)/(outside diameter 
of pipe) 


tip rotation 


xu/L 


ee 


he b+ %0h2+%0 
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SECTION II 


THEORETICAL ANALYSIS AND TECHNIQUE FOR INVESTIGATION 


2-1 Review of Flexibility Analysis for the Ell-Shaped Pipe 

Before presenting the investigation of this ell-shaped piping 
configuration in detail, it is felt that certain principles re- 
garding the analysis merit some review. 

If the tip loaded 
cantilever beam shown in M 
Figure 1 is analyzed by 
writing a moment equation 


about point A the following 


equation results: oS L ab 


Figure 1 





oe ETy = M + Q(d-y)+ P(L-%) 


Solving this differential equation and applying boundary 


conditions yields 


(2) y= MiZul- cost secu + PCI rN ub sin §+ fic cos) anu] : pest 


Sib Ue 
for boundary conditions a. y(0) = 0 
b. y(L) = A 
ce. y'(0) = X 
d. y'(L) = P 


See Table I for nomenclature. 


In addition 


G) A = (Ma, + PL(o, + 70,)](¢/er) 


14 


and 


(4) De [ Ma, + PL(%,+ ) oly) ( L /et) 


The alpha coefficients appearing in the above and subsequent 
equations are defined in Table II with alternate forms to facilitate 
computation, 

In terms of all four secondary effects, 6is given by equation (5). 

2] 
6 §& = eL- 2vplp — (Q-Acb)L ~ PY 
E()-p*) AE 

The first term on the right hand side of equation (5) represents 
thermal expansion, the second represents axial contraction resulting 
from hoop tension caused by internal pressure, the third represents 
axial contraction due to axial compression in the metal, and the 
last represents the effective axial contraction resulting from the 
difference between arc and chord. The pressure terms may be combined 


to yield the equation, 


(6) & = (e+pH)L - QL/AE - A 


where y is given by L 


= | \\e 2 
(7) \ = z \ ((yP-X*) ax 
° 
when higher order terms are neglected (see Appendix G). If alpha 


coefficients are substituted, 


(6) A= (MPdy + 2.0MPLaly + P*ZX,.)( U/er) 


15 


Other pertinent relationships include the shear distribution 


eaecor- 


(9) G = 6, + He PY (| + p?)* 


where 


ao) 


i 


(74 140 48V )fo (140) 


and 


ay fAe= (10 #20V4 8*)/a(i +0)" 


To apply the foregoing analysis to an ell-shaped pipe, the 
corresponding quantities are defined as they appear in that 
configuration. The basic structure is that shown in Figure 2. If 
the configuration is subjected to uniform heating, thermal expansion 


creates a situation resembling that of Figure 3. 


34 
jt 


a x-bF 


Figure 2 Figure 3 





The resulting free body diagramsfor the two legs a and b is depicted 


in Figure 4. 


John E. Brock, "Shear Distribution in Piping," Heating Piping and 
Air Conditioning, Vol 35, January 1963, pp. 141-143. 
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Figure 4 


Referring to the horizontal leg as part 1 and to the vertical leg as 
part 2, Table III list the quantities corresponding to the general 


case for each part. 


TABLE III 
ITEM PART 1 PART 2 
L a b 
P -F F, 


“RQ CnC = o 
> 
K 


Geometrical continuity requires equating corresponding expressions 


for X,Y, and @; i.e.; é, > ive 
bv = Ao 
P= Po 


17 


or 


(iz) @a~ Faye - (A+ bak = [-NG%, + Rb (A. + 44;)} béT 
(13) @b- E, Hag. - | +pbH =[-Ne + Fa(ot, 49 %)| ky 
14) Naty Fab (o,+ 94) = - Nba, + F b°(&, + 5%) 


A superior bar indicates that the quantity has been evaluated 
using the parameters of part 2 and no superior bar indicates evaluations 
with the parameters of part l. 


Solving equation (14) for the moment, N, yields 
> (= * 
(is) Ne Far (deny \ + Fab* @ +7 &, 
AAs + bX, 


where (15A) ‘ey = a” (+o) fa. ol + bos) 
sala sb (5, +5.) fava, +b) 


Now define; 


a6 d= -%, 8 EI 
(17) ie =[b Ol, + al, ) is at —o 2b ft 


wo | ted i194) pL -o, 90°} /e7 
(19) = A, Pir 


and the system of equations become 
di 4 F (e+pk) a -h 

(20) —y 
ai dn} [| AJ™ L(espeb —A 
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In order to solve for the inknowns Fo and Po» an iterative 
procedure is strongly suggested by the nature of the equations. 

By assuming initial values of Fy and F, »which can be any 
value, one can calculate finite values for all the coefficients 
85 j and the quantities on the right of equation 20). This equation 
is solved for new values of Fy and i) which in turn provide improved 
values for the next iteration. This procedure of iteration is 
continued until satisfactory convergence is reached. Once convergence 
is obtained such quantities as moment, 9, and other desired values 
are obtained using the converged values of Fy and F,: 
2-2 Dimensionless Parameters 

The large number of variables involved in this problem could 
produce millions of different problems. In order to pursue the 
investigation as thoroughly as possible and at the same time retain 
some compactness in the results, four dimensionless parameters were 
used, These are: 

(i hee Ratio a/b 

(2) (Radius of Gyration)/(Total Length) 

(3) Pressure/Modulus of Elasticity 

(4) Thermal Strain 

To cover reasonable ranges of all four parameters evidently 
requires solution of a large number of problems each of which is 
itself rather demanding. Accordingly, use of a large digital computer 
is definitely required. 

Each of the four secondary effects can either be included or 


excluded in the calculation, This produces is or 16 different possible 


combinations of secondary effects. 


lie 


2-3 Requirements Imposed on Digital Computer Program 

A computer program for studying this problem has certain basic 

requirements. 

A, First of all,it must be capable of accepting certain input 
data, which includes a basic physical description of the 
problem as well as control variables that instruct the 
program in how to treat these data. 

B. From the input information it must recognize the parameter 
to be investigated and must load the various values of the 
parameter over the desired range. 

C. For each parameter value in the given range it must: 

1. Solve for the F's by the iteration procedure and then 
determine other desired values. 

2. Accomplish C.1l. for all different combinations of the 
four secondary effects. 

D. It must increment the parameter value and repeat step C 
until results for all values in the range have been found. 

E. It must provide results in either tabular or graphical form. 

The following section describes the computer program in more 


detail and shows how the above listed demands are met. 
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SECTION IIT 


DESCRIPTION OF COMPUTER PROGRAM 


3-1 General Remarks 

As mentioned in Section II and suggested in reference 2, a 
computational procedure is evident. Assume initial values of 
Fy and Fy (which may be taken to be zero), evaluate the coefficients 
Bj and the terms (and A and upon @lving the system of equations (20), 
obtain improved values of Fo and F,> continuing until satisfactory 
convergence is reached. This idea along with the other basic require- 
ments (A-E of Section II) composes the foundation upon which the 
computer program is constructed, 

The program was first written in Fortran 60 computer language 
for use on the Control Data Corporation Model 1604 digital computer 
located at the Naval Postgraduate School, Monterey, California. In 
the midst of the investigation reported in this paper, the CDC 
computer was replaced at the Postgraduate School by an International 
Business Machine Model 360/67. The present version of the program 
is written in Fortran IV language for use on the latter machine. 

All computations in the program are executed using the Double 
Precision Capability of the system. This feature improves the 
accuracy of the results but unfortunately requires more time than 
the single precision calculation, 

An effort was made to keep the program as general and as flexible 
as possible. All of the information needed to specify a given problem 
is read directly from data cards. The information on the data cards 
defines the physical condition to be investigated and gives the program 


additional control information concerning the parameter to be used 


2: 


in evaluation, range of the parameter, combinations of the secondary 
effects to be investigated, convergence criteria, etc. 

It allows the user four different parameters with which to 
evaluate the problem. It allows for all sixteen possible combinations 
of the four secondary effects and computes all sixteen of these for 
each value of the parameter. 

It has optional output. In the data cards the user specifies 
the desired type of output: tabular, graphic, or both. 

To describe the computer program more completely requires 
explaining in more detail the functions of its seven major parts. 

A. Input 

B. Loading of Parameter Values 

C. Consideration of All 16 Different Combinations 

D, Computation by Iteration 

E. Incrementing Parameter to Next Value in Range 

F. Output 

G. Subroutines 
3-2 Input 

The term input is used here in referring to the information that 
is supplied to the program via the data cards. In all, there are 
fifty-five different items read in as data. These items completely 
define the physical situation under consideration. Further, it informs 
the program of the dimensionless parameter to use in the investigation, 
the range of the parameter, conditions to investigate, and so on. 

Each of the variables read in as data is listed below along with a 
description of its function in the program. The order of this listing 
is the order in which they appear on the data cards. The names of the 


variables are the names which are used as input. The same physical 


Le 


quantity may be called by other names later in the program and it 
should be remembered that the variable name mentioned below is 


applicable only in the input section of the program. 


R length of the shorter leg, a, (inches) 

Ss length of the longer leg, b, (inches) 

ZIN pipe inside diameter (inches) 

ZOUT pipe outside diameter (inches) 

DTEMP temperature difference (degrees fahrenheit) = 
el ; ares ai” 

Q internal pressure (psi) 

E modulus of elasticity of pipe material (psi) 

COEFI coefficient of linear expansion of pipe material 
(in/in - °F) 

KPARL Integer variable which can assume the integer values 


1 to 4 inclusive, that designates the parameter 


to be used in the investigation as follows: 


KPARL = 1 = Leg Ratio 
KPARL = 2 = (Radius of Gyration)/(Total Length) 
KPAR1 = 3 = Pressure/Modulus of Elasticity 
KPARL = 4 = Unit Thermal Strain 
L Integer variable that specifies the size of the major 


parameter increment (see Section 3-3). This can be 
any value but for graphic output should not require 


more than thirty stored values of the parameter. 
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Integer variable that designates the initial value 
of N, itself an integer variable which is used 
as a subscript for the parameter values. This 
(the value of NN) is useful if the program ever 
exceeds the time allotted for one run and it is 
terminated by the operator, or if the given size 
of increments do not yield convergence over the 
desired range. By giving NN the value of the last 
value of N on the terminated run, the user is able 
to modify certain other control variables such as 
L, LL1, FTE's, etc., for the next computer run and 
begin computation at the point where the last run 
ended, regardless of the reason for termination. 

For the first run of a given parameter or problem 
the value of NN is always 1. 

Integer variable that designates the final value of the 
subscript N in the parameter range. 

Integer variable which designates the initial value of 
MT. The significance of MT is explained in more 
detail in the section defining loading. MB can take 
on any value from 1 to 16. 

Integer variable that indicates the final value of MT. 
MC = MB always and 1S&MC<16. Both of the above 
conditions must be satisfied. The program is designed 
in such a way that if MT does not equal its final 
value, MC, then MT will be incremented by 1, or the 
next combination of secondary effects will be 


calculated using the same parameter value. What 
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MIT 


KKK 


In 


MD 


MD 


MD 


ibe 


In 


this in effect does is have the computer make all the 
calculations for all the combinations between and 
including MT = MB and MIT = MC. For example, if 

MB = 3 and MC = 7 for a given set of input values, 

MT would take on values 3, 4, 5, 6, and 7 for each 
value of the parameter. Consequently, the output 
would reveal results for the combination correspond- 


ing to the 5 different values of MT. 


teger variable that designates the type of output 
desired, 
= MD #*3 
= 1 Tabular output only 
= 2 Graphic output only 
= 3 Both tabular and graphic output 
teger variable that indicates the limit to the 


number of iterations the program will perform. If 
this number is reached before the convergence 
criterion is met, the program will terminate 
operation on its current set of data, give the 
necessary output, then proceed to the next set 

of data. 

teger variable that assigns a value to the constant 
IA. IA is used in a relationship that determines 
the improved values of F's to be used in the next 
iteration on the basis of the last assumed values 


of the F's. The relationship involved is: 


Fis 4x Fl(New)+ LAx Fi (ow) 
440A 


Zo 


LL1 


EE2 


EXSCAL 


YSCALE 


IXUP 


LYRIGH 


FTIE11L=-FIE 162 


and a similar relationship exists for the improved 
F2. It was found that generally as the value of N 
increases IA should increase also, thereby making 
the improved value closer to the old value than 
the new value. 

Integer variable that specifies the size of the sub- 
increments of the parameter. (See section concerning 
loading). 

This value sets the convergence criteria, After each 


iteration the new values of the F's are compared to 


the old ones using the equation 





If this is satisfied the values are considered converged. 
This quantity is the argument used in the DRAW subroutine 
(see Appendix A-4) that assigns the scale to the X-axis 
or abscissa axis when graphic output is desired, 
Same a EXSCAL above except it governs scale on Y-axis. 
Another argument for the DRAW subroutine that specifies 
how many inches up the graph one wants the X-axis 
located, 
Specifies how many inches from the left side of the graph 
one wants the Y-axis. 
The 32 different values are the last converged values 
of Fl and F2 for a particular combination of secondary 
effects. The final digit distinguishes between Fl 
and F2. The other numerical digits run from 1 through 


16 and correspond to values of MI. These FTE values 
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are retained by the computer and uSed as the first 
approximation for the F's in the corresponding 
combination when the parameter has been incremented. 
In most investigations the values assigned by the 
data card are all zero. These values, however, are 
just the initial values assigned to the FTE's. Once 
computation begins the computer takes over the task 
of assigning these values. Every time convergence is 
reached the values of two FTE's change and as 
convergence is reached for all combinations and one 
parameter value, all 32 change their walues. Use of the 
variables aids the computation process in two ways. 
The nature of the governing equations causes difficulty 
with convergence under some conditions, Therefore, 
if the first assumed values of the F's are near the 
correct values a better chance for convergence exists. 
If, on the other hand, the guess is arbitrary, or not 
close to the converged value, the system of equations 
may not converge. If the last converged values are 
assumed as the first approximations with the incremented 
parameter, they should be good estimates depending on 
the size of the parameter increment. The smaller the 
parameter increment the better the estimate should be. 
Another utilization of the 32 variables allows 
the user the capability of picking up computation at 
any point in a parameter range. Suppose, due to excess 
time the run was terminated before the entire range 


was completed, Then, on the subsequent run, the user 


Zi 


would be able to start at the same place in the parameter 
range that the previous run ended, This provides 
continuity without duplication and the ability to read 
reasonable first approximations to the F's in as data 


regardless of where in the range one desires to continue. 


3-3 Loading of Parameter Values 

The integer variable KPARI1, read from the data card, instructs 
the program which dimensionless parameter will be used to investigate 
the Seer e Cure: For a given set of data (one set of six data cards) 
only one parameter can be investigated. Each of the four values that 
KPAR1L may take on indicates a different loading scheme for the six 
ected variables. Once the program knows the value of KPARI it 
proceeds to the proper loading scheme for that parameter. 


The six sSubscripted variables involved are: 


AA(N) leg a 

BB(N) leg b 

XDIN(N) pipe inside diameter 
XDOUT (N) pipe outside diameter 
XP(N) internal pressure 
ALIT(N) unit thermal strain 


Here again the variable names associated with the physical 
quantities apply only in the loading scheme. The same physical quantity 
may have another variable name at another place in the program. See 
Appendix A for a complete listing of computer program nomenclature. 

The subscript N takes on integer values from 1 to MA. For each 
value of N, the six subscripted variables have values stored in the 


computer memory. In each of the different loading schemes only one 
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of the six subscripted variables takes on different values for each 
different value of N. The other five have fixed values. 

Associated with N in the loading procedure are two other integer 
variables, MA and L. The values of these two quantities determine 
the range of the parameter and the size of the main parameter increment. 

An important feature of the loading process is that the parameter 
values are loaded in such a way as to produce the least critical loading 
conditions on the structure for the parameter value associated with N = 1 to 
the most critical loading condition when N = MA, This essentially, 
means that at the least critical end of the range the initial assumption 
of the arbitrary values of the F's will probably yield convergence. How- 
ever, as the parameter proceeds to the more critical end of the range, 
the facility utilizing the FTE's becomes more and more necessary in 
achieving convergence. 
3-4 Combination of Secondary Effects 

The capability that allows the computer program to solve the system 
of equations for all sixteen different combinations of secondary effects 
for each parameter value is afforded through the subroutine FACTOR. 
To supply the subroutine with information it needs, the integer variable 
MIT is used, MT takes on values 1 ~ 16 for each parameter value in the 
range. Each of the 16 values represents a different combination of 
secondary effects, When the subroutine is called with MT specified, 
it responds by supplying the main program with four factors and a 
particular designator. These four factors are multiplied by the terms 
in the system of equations that represent the four different secondary 
effects, Each of the factors can have only two different values: 1.0 
and 0.0. It is therefore evident that if the combination requires the 


suppression of an effect in the calculation, the governing factor is 


Zo 


zero; on the other hand if the effect is to be included the 
governing factor is unity. 

The main program is so designed that all of the sixteen different 
combinations do not have to be calculated each time, say in a situation 
where the user is only interested in one or two combinations. If such 
is the case the integer variables MB and MC are manipulated to provide 
the combinations desired. MB sets the initial value of MT. After 
each combination is completely solved for all desired quantities the 
main program tests the value of MT against that of MC. MC represents 
the final value of MT to be considered for each parameter value and 
in the test if MT is less than MC its value is incremented by 1 and 
that combination is investigated, This process continues until 
MT = MC and after the MT = MC combination is solved the program increments 
the parameter. The limitation on this capability lies in the fact that 
combinations represented by values of MT between MI = MB and MT = MC 
will be calculated even if some of these intermediate cases are not 
specifically wanted. 

The designator referred to earlier is realized as a label in the 
printout of the results. It is printed in a certain column of the 
results and is a four character identifier that tells the value of 
each factor in that particular combination (see the explanation of this 
subroutine for specific meanings for these designators in Section 3.8). 
3-5 Computation by Iteration 

The iterative procedure, as mentioned earlier, is that of assuming 
initial values of Fl and F2 (which may be taken to be zero), evaluating 
the coefficients Bi; and the terms A and > and upon solving the system 
of equations (20), obtaining improved values for the F's. This is 


continued until satisfactory convergence is reached. Through the use 
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of the values of the FTE's an attempt is made to minimize the number 
of iterations required. This saves computer time as well as increases 
the probability of convergence. 

After each iteration the newly calculated values of the F's are 


tested against the old ones using the following convergence criterion: 


| Fi(New) — F L(ovo)| + F 2 (NEW) - F2 (o.0)| 


= EPC 
Fi(vew))” + (F2(NeEw))* 


If this criterion is not satisfied, the program selects improved 

values for the next iteration and sends these newly selected assumptions 
of Fl and F2 through the system of equations again. The process for 

the selection of the new F's is defined by the following two relation- 
ships: 


FINEXT= 4x FIN+IA*FL. F2NEXT= 44 F2N+IAX F2 
4+IA ? 4+1A 


where FINEXT and F2NEXT are the improved values of Fl and F2 respectively 
to be used in the next iteration. FIN and F2N are the new values of 
Ey and EF, just computed and Fl and F2 in the equations are the values 


of EY and E, assumed at the beginning of the iteration. The factor 
IA is an integer whose value is set by the integer variable, KKK, 

read in from the data card. It is obvious that the improved value can 
be closer to the older value or the newer value at the discretion of 

the user. The nature of the equations is such that when the parameter 
value produces less critical loading conditions, convergence is achieved 
more rapidly if the improved value is taken to be closer the new values. 


Conversely, when the parameter is in a more critical region, a selec- 


tion of an improved value closer to the old value produces convergence 


oi 


more rapidly. The term "critical region" 


refers to that range of the 
parameter where the system of equations behave in a highly nonlinear 
fashion, 

There is no set rule for this selection of improved values and 
the best instruction is gained through experience with the program. 
Using values of IA between 0 and 10 is normally adequate for most 
Situations. 

If the convergence criterion is satisfied, the latest calculated 
values of the F's are accepted as the correct ones for the given value 
of the parameter and the given combination of secondary effects. 

Using these values of the F's other necessary quantities are calculated. 
The program then goes to the next combination of secondary effects 

and begins iteration for the solution of the F's for that situation. 

The first approximations for the F's are the converged values found 

for that combination with the previous parameter value. When the 
program has solved the equations for all combinations desired, it 
increments the parameter and the process starts all over again. 

3-6 Incrementing the Parameter 

Two different parameter increments are utilized in the program. 
They are referred to hereafter as APAR and Spar. The two can be 
described as follows: 

If PAR(N) is the parameter for a particular value of the subscript, 
then PAR(N + 1) = PAR(N) + QPAR. Therefore, AVPAR is the increment 
whose size is the difference between two consecutive stored parameter 
values. The size of QPAR is determined by appropriate input data 


and by the control variable L (see description of L in Section 3-2). 
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The sub-increment, Opar is used to give the program the capability 
of making the parameter as small as deemed necessary in order to achieve 
convergence. The basic theme of the investigation of any of the four 
dimensionless parameters is to proceed from that end of the range 
that provides the least critical load toward the end which produces 
the most critical condition in the range. Each time the parameter 
is incremented, the last converged values of the F's are used as 
the first approximation. By making the increments smaller, the first 
approximation is closer than if the parameter had a larger increment. 
Therefore, it is conceivable that by making these increments small 
enough, convergence can always be realized, It is for this reason 
that the program uses increments of two different sizes. 

For the graphic output capability, it is necessary to store many 
values of the variables that are used in the graphs. The storage 
capacity of the computer thus establishes a lower limit for the larger 
increment. In order to further decrease the size of the increment the 
smaller increment, Spar, is used. Except when they coincide with 
APAR values, parameter values resulting from this smaller increment 
and the corresponding calculated values are not stored, but used only 
to obtain new FTE's for the subsequent parameter value. Another way 
of expressing the process is to say that after all necessary calculations 
are made for one parameter, the parameter value is incremented by Spar. 
Using the last converged values of the F's as first approximations, 
converged values for the new parameter value are sought. Thus the 
incrementation continues over the parameter range. If the first 
parameter value considered was a stored value, then after LL1 sub- 
increments the next stored parameter value will be reached. LLl is 


the integer variable read from the data card. 
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Figure 5 aids in understanding the use of the two increments. 
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Figure 5 


From Figure 5 the following relationships are apparent: 

LL1 x Apar = APAR PAR(N + 1) = PAR(N) + QPAR 
Also, for the parameter values corresponding to the solid vertical 
lines, necessary calculated values are stored for graphic or tabular 
output. For parameter values corresponding to the broken vertical 
lines only converged values of the F's are calculated for use as 
the first approximation when the parameter is incremented. 

To summarize, each time the parameter is incremented it is 
increased by the amount Spar. Every LL1 times, all calculated values 
necessary for output are stored. The size of APAR is somewhat 
restricted by the storage capacity of the computer but the size of 
6 par has no restrictions, It must be remembered, however, that by 
decreasing the size of } par , the computing time for a given range of 
the parameter is increased accordingly. Therefore, in the interest 
of time on the computer, the value of dpar should be as large as possible, 
commensurate with providing the necessary convergence. 

3-7 Output 
The user has some choice concerning the type of output. The 


results can be represented in tabular form, graphical form, or both. 
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To specify the type of output desired for a given set of data,assign 
the proper value to the input integer variable MD. MD can take 


on values 1 to 3 with the following significance attached to each: 


MD = 1 Tabular Output only 
MD = 2 Graphic Output only 
MD = 3 Both Graphic and Tabular Output. 


If MD is set equal to 1, ten columns of calculated values will 
be printed out for all of the desired combinations of secondary 
effects. The values shown are: the value of the subscript (N), 
the value of the integer MI, the combination designator, moment, 
moment ratio, the converged values of F, and Fy» the value of the 
parameter, and the number of iterations required for convergence 
in each case. In addition at the top of each page of output are 
listed the values of the raw data exactly as it is read in from the 
data card. These values only serve as a convenience in showing what 
was read into the program and what the individual values were prior 
to the parameter loading. Other values that are printed include the 
exact moment, the parameter that was investigated, the elapsed time 
in seconds since the last stored parameter value with its associated 
values were calculated, and the critical loads for each parameter 
except the leg ratio. Actually the output sheet is very straightforward 
if one follows the column labels. 
Some of the aoe mentioned terms deserve a little more explanation. 
Combination Designator: The four character designator supplied 
to the main program by the subroutine Factor. The four characters are 
either zero or one and they represent the value of a particular factor 


for the given combination of secondary effects. For instance, if the 
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designator is 1001, the factor that is multiplied by the term represent- 
ing the secondary effect, arc and chord difference, is 1. The factor 
that is multiplied by the secondary effect term for shear is 0. The 
factor multiplied by the term representing axial deformation is O and 
the factor multiplying the beam column term is 1. The order of the 
factors always appear in the same order. Arc and chord, shear, axial 
deformation, beam column action appear in that order from left to right. 

Exact Moment: When the term exact is used referring to a calculated 
quantity in the system of equations it means that the calculation was 
made with all four of the secondary effects allowed and none of the 
effects suppressed; for such a calculation, the combination designator 
1S EL. 

Moment Ratio: Ratio of moment calculated with given combination 
designator to exact moment. 

If MD is set equal to 2, only graphic output results. For graphic 
output certain limitations exist: 

a. The ordinate is always the moment ratio. 

b. The abscissa is always the parameter under investigation. 

c. MA or the number of stored values for each subscripted 
variable should not exceed 30. 

One graph is obtained for each set of data investigated. On each 
graph six sets of values are plotted. These six sets correspond to the 
first six values of MT. See Subroutine Factor in Section 3.8 for the 
combinations of secondary effects that correspond to the six values. 

The combinations included in the six are with all effects allowed (exact 
solution), all effects suppressed, four combination where only 1 of 
the four effects is allowed and the other three are suppressed. There- 


fore, if the user specifies graphic output only, the input integer 


B16; 


variables MB and MC should equal to 1 and 6 respectively. Calculations 
for all other combinations would be redundant since there is no 

tabular output and the first six combinations are all that will be 
represented on the graph. If on the other hand, tabular output is 
selected or the user desires both types of output on the same run 

then the above restriction is not applicable. 

Before one employs the graph plot form of output it is necessary 
to have some knowledge of the values that are to be plotted. The 
data card contains four of the arguments to the DRAW subroutine. 

For a meaningful graphical output a reasonable scale has to be 

selected and other parameters have to be chosen. The recommendation 
therefore is to use the tabular output first and from it plan the 
graphical output so as to describe the results most usefully. [In 
addition the user should be familiar with the information relating to 
the DRAW subroutine contained in the next section and in the subroutine 
system library write-up reproduced in Appendix A-4, 

3-8 Subroutines 

There are three subroutines used in the program. The first two 
subroutines, designated ALPHA and FACTOR were programmed by the author 
and the third, designated DRAW, is from the subroutine library at 
the computer facility of the Naval Postgraduate School, Monterey, 
California. The complete library write-up forDRAW and its use appears 
in Appendix A-4, 

ALPHA: This subroutine is used by the main program to compute 
the values of the alpha coefficients as they are shown in Table IT, 
page 13. It is called twice during each iteration and supplies the 
program with values of the As orQ$ corresponding to the information 


that is given to it, namely, &J, ) ratideW eI or ws?) , and EPI, 
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depending on whether the main program wants the valuesof o's 
or OA’, . 

FACTOR: This subroutine is used by the main program to set 
the combination of secondary effects to be considered. As previously 
mentioned, it works along with the integer variable MT which takes 
on values 1 through 16. Each of the values designates a certain 
combination of the secondary effects. Upon receiving this number 
designation, the subroutine supplies the program with four factors 
each of which has the value of either 1.0 or 0.0 depending on the 
combination. These factors are multiplied by the quantities that 
introduce the secondary effects into the system of equations. If 
the factor is 1.0 that particular secondary effect is included in 
the calculations; if it is 0.0 the secondary effect is suppressed in 
the calculation. Tables IV and V on the following two pages give more 
information concerning these combinations and their designations 
that is helpful in understanding the particular system for varying 
the combinations. 

DRAW: This subroutine is the one through which the program 
affordsthe user the graphic output option. The complete write-up 
of this subroutine just as it appears in the Postgraduate School's 
subroutine library is reproduced in Appendix A-4. It is felt that 
no further explanation is necessary except for the four arguments 
that are read in as data in the main program. Referring a the 
write-up in Appendix A-4, the arguments h. through k. are supplied 
to the calling statements by the last four values on the second 
data card. This allows user control of the scales desired and the 


axis location. 
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TABLE IV 


FACTOR NAME SECONDARY EFFECT THAT IT CONTROLS 
FAC1 Difference in Arc and Chord 
FAC2 Shear 
FAC3 Axial Deformation 
FAC4 Beam Column Action 
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A means effect iS allowed 
B means effect is suppressed 
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SECTION IV 


RESULTS OF INVESTIGATION 


4-1 Errors in Paper No. 63--APMW-21 

Prior to developing the computer program to utilize the relation- 
ships of Section II, certain review and checking was undertaken by 
the author. This review was undertaken principally to gain a working 
familiarity with the equations and nomenclature. During the course 
of this review and later in developing the program a few errors 
were detected in the original paper. 

Most of these mistakes were algebraic in nature. One was in 
the area of the critical buckling loads for beam column considerations. 
The mistakes are: 

1. Equation (8), the factor "2.0" in the second term of the 
numerator did not appear in the same equation in Brock's paper. 

There were three mistakes in the formulas for the Of coefficients 
of Table II. 

2. & 3. In the formula for oO ; the &'s in the first two terms 
of the series should be ol, 's. In the reference paper they are of, 5 : 


4, The formula for @ 


10? the last term should be mF ot as 


opposed to "ot al '" in the paper. This was apparently a typographical 
error. 
>. The last incorrect item found was in the statement relating 
iy 
to the critical buckling loads. The paper indicates that &f= Wa. 
defines the critical buckling load for the ell-shaped pipe structure. 
This is true only for the equal leg case but is not true in general. 
All of the results reported in this paper reflect the corrections 


as stated above. 
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4-2 Numerical Calculations 

With the use of the digital computer program described in 
Section III, calculations were made for several sizes of pipe 
and for various values of temperature increase, as well as for 
different values of the other inputs. This was done not only to 
test the integrity of the computer program, but also to obtain 
sufficient data to permit a search for trends and useful modes 
of presentation. It turns out that the leg ratio, a/b, is the 
most significant parameter, as might have been expected, and 
presentations are confined to displaying the moment ratio as a 
function of leg ratio for various combinations of secondary 


effects all curves shown being for the following specific case. 


Length of leg a varied 

Length of leg b 360 inches 

pipe inside diameter 10.020 inches 

pipe outside diameter 10.750 inches 
Temperature difference 530 degrees fahrenheit 


(70-600 degrees F) 


pipe internal pressure 400 psi 

pipe modulus of Elasticity De 10° psi 

pipe coefficient of linear 7325 x iin in/in deg F 
expansion 


4-3 Final Results and Conclusions 

The results are shown in the form of eight curve sheets, 
Figures 6 through 13, which should be considered in pairs, the first 
of each pair showing the effect of suppressing (i.e. failing to 


consider) only a single secondary effect and the second showing the 


4.2 


effect of suppressing the other secondary effects and considering 
only this one. For comparison purposes, each curve sheet also 
shows the effect of neglecting all four secondary effects. 

These curves are based only upon what has been called the 
‘moment ratio", i.e., the ratio of the moment at the intersection 
of the two legs as calculated by an inferior procedure to that 
calculated by the best available procedure, namely that of including 
all four secondary effects. It is clear that substantially the same 
sort of result would have been obtained by considering a similar 
moment ratio for either of the two built-in ends or a similar ratio 
of reaction force. 

These curves may be regarded as speaking for themselves. The 
most significant conclusion that may be made is that neglecting 
secondary effects, indeed can have a profound influence on the 
accuracy of piping flexibility calculations, leading to gross errors 
in the evaluation of stresses in cases of leg ratios which, though 
extreme, are likely to be encountered in piping practice. 

Since the secondary effects are quite subtle to evaluate, this 
conclusion could not have been forthcoming from a less sophisticated 
and demanding approach. It is believed to be an entirely new cone 


clusion and it is certainly a significant one. 
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APPENDIX A 


COMPUTER PROGRAM NOMENCLATURE 


A-1 MAIN PROGRAM 

NOTE: Where possible the variable names of the computer program 
are identified by the notation of Tables I and II. Also, in some 
cases the same physical quantity may be represented in the program 
by more than one variable name. In such cases the extreme right 
hand column will indicate the general area of the program where the 


physical quantity takes on the indicated variable name. 


VARTABLE NAME TABLE I & II NOTATION COMMENT 

A a Computation Section 
AA (N) a Loading Section 
A1-A10 OX = Aro 

AB1-AB10 X, — Ayo 

ALIT(N) e Loading Section 
ALPHA Subroutine ALPHA 

AREA A 

B b Computation Section 
BB(N) b Loading Section 


Bl (erpHa - A 
B2 (e+pH)b- A 
vent La (et, + 94) ]/(a.04, + bats) 
BET CO (a,+)%,)]/(a0,+ bB;) 


COEFT Material coefficient of thermal expansion 
DELT1-DELT4 Size of 8 par for each parameter 1-4 respectively 
ie JuxkQe- Jiux dS. 

DIN pipe inside diameter Computation Section 

DOUT pipe outside diameter Computation Section 
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COMPUTER PROGRAM NOMENCLATURE (cont.) 


VARIABLE NAME TABLE i & IT NOTATION COMMENT 

DTEMP Difference of Temperature 

E E 

ELIT e Computation Section 
EP1,EP2 G, E> 

ETA n 

ETAB y) 

EXMOM (N) Exact moment for a particular parameter variable, i.e. 


with all secondary effects allowed. 
EXSCAL Applicable to graphic output: X scale in units per inch 
FAC 1-FAC4 Four factors used for either suppressing or allowing 


their corresponding secondary effects. 


FACTOR Subroutine FACTOR 
Biee2 Assumed values of the forces FF, at the beginning 
of an iteration. 

FICRIT,F2CRIT Euder Buckling load for Fi and Fo. 

FIN, F2N Newly calculated values of Fy and F, from a particular 
iteration, 

FTE11-FTE162 See section on input Data 

G G Computation Section 

85; coefficients defined in equations (16)-(19) 

H H 

II Counter use to count the number of iterations required 
for convergence, 

ITA Factor to aid in determining improved values of F's 

LABEL Argument for Subroutine DRAW 

ITITLE Argument for Subroutine DRAW 
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VARIABLE NAME 


IXUP 


LTYRIGH 


KHAT 


KKK 


KPARL 


LABEL 


LL1 


MC 


MIT 


MT 


NN 


PAR2 


COMPUTER PROGRAM NOMENCLATURE (cont.) 


TABLE I & It NOTATION COMMENT 
Argument for Subroutine DRAW 
Argument for Subroutine DRAW 
Counter used to count the number of subincrements between 
the larger main increment. 
Secondary Effect combination designator 
Input variable that assigns value to IA 
Input variable that designates parameter to be investigated 
Input variable that aids in determining QPAR. 
Argument for Subroutine DRAW 
Input variable that sets size of subincrements, épar. 
Input variable that designates final value of subscript,N. 
Initial value of MT 
Final value of MT 
Input variable that designates type of output 
Input variable that specifies the maximum number of iterations 
that will be performed before terminating operation on a 
particular set of data. 
Integer variable that specifies the combination of secondary 
effects to be considered. 
Subscript, used to identify calculated values of desired 
quantities with the corresponding parameter value. 
Initial value of the subscript N 
Pp Computation Section 
A designation for the second parameter, (Radius of 
Gyration)/(Total Length), in the computational part 


of the program. 


VARIABLE NAME 


PIE 
PIESQ 


POLIS 


Q 


WBAR 
XDIN (N) 
XDOUT (N) 
Salen 
XLAB 
XMU 1 
XMU2 

XN 
XP(N) 
YSCALE 
nat 

ZIN 
ZLEGR(N) 
ZLEN 


ZMOM 


COMPUTER PROGRAM NOMENCLATURE (cont.) 


TABLE I & II NOTATION 
TT 
T° 


P 


a 


Radius of Gyration 


es pe 2 wm 


pipe inside diameter 


pipe outside diameter 


/\ 


Va 
Fz 


COMMENT 


Computation Section 
Input Data Section 


Input Data Section 


Input Data Section 


Loading Section 


Loading Section 


Moment in Computation Section 


P 


Applicable to graphic output: 


I 
pipe inside diameter 
First parameter, a/b 
Total Length, a +b 


moment 


56 


Loading Section 


Y scale in units per inch. 


Input Data Section 


COMPUTER PROGRAM NOMENCLATURE (cont.) 


VARIABLE NAME TABLE I & II NOTATION COMMENT 

ZOOM Initial thermal strain 

ZOUT pipe outside diameter 

ZPAR2 (N) Second parameter, (Radius of Gyration)/(Total Length) 
ZPAR3 (N) Third parameter, Pressure/Modulus of Elasticity 


ZRAT1(N)-ZRAT16(N) Ratios of the Moments of the sixteen combinations 
divided by the exact moment for that combination. 


ZRATIO Same as the ZRAT()'S except this value is not stored. 


A-2 SUBROUTINE ALPHA 


VARIABLE NAME § TABLE I & II NOTATION COMMENT 
. yor] 

EPL €, 

T1-TLO OHy—- Aj OF K, — Hp 

TX lwl er lw] 

(-w)*or (-ax)* 

Z (us)* or (w)* 


A-3 SUBROUTINE FACTOR 

VARIABLE NAME TABLE I & If NOTATION COMMENT 

JT Integer variable, MT 

LLAT Secondary effect combination designator 

R1-R4 Four factors used for either suppressing or allowing 


their corresponding secondary effects. 
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A-4 SUBROUTINE DRAW-DESCRIPTION FROM SUBROUTINE LIBRARY 


A, 


IDENTIFICATION: 
TITLE; General Graph Output Subroutine 


SUBROUTINE NAME: DRAW(J7-NPS-DRAW for CDC 1604 revised for IBM 360) 


DATE: July, ivory 
PROGRAMMER : Revised by Patricia Johnson 
PURPOSE : 


This subroutine generates data for plotting graphs on a CalComp 
plotter. Provision is made for curve drawing and point plotting, 
automatic scaling, graph titling and axis annotating. An attempt 
was made to provide a considerable amount of flexibility, at the 
expense, necessarily, of a relatively large number of arguments 
and a rather high memory requirement. 
USAGE : 
1. Definitions: 
In what follows the word "graph" will be taken to mean one 
piece or frame of graph paper on which there may be plotted 
one or more curves and/or sets of points. A "curve" will 
mean a continuous line generated by the sequence of straight 
lines joining successive points of the set defining the 
curve. A "point plot'' will describe the representation of 
a succession of points by means of symbols (such as a cross) 
on the graph. The points are not connected ina point plot. 
2. Calling Arguments: 
All necessary information is transferred to DRAW through the 
calling arguments. The call statement is: CALL DRAW (NUMPTS, 
X,Y ,MODCUR , ITYPE , LABEL, ITITLE , EXSCAL, YSCALE , IXUP, ITYRIGH, 


MODXAX ,MODYAX , IWIDE , IHIGH, IGRID, LAST) 
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It is important to realize that one and only one curve or 


set of points is plotted each time DRAW is called. However, 


it is possible to call DRAW repeatedly if several curves 


and/or sets of points are wanted on one graph. 


The calling arguments are as follows: 


a. 


d. 


NUMPTS : 


MODCUR: 


The number of points defining a curve (28 NUMPTSS 900), 
of the number of points to be point plotted (2% NUMPTS$30). 
The array of X-ordinates, X, must be dimensioned at 

least equal to NUMPTS in the calling program. 

The array of Y-ordinates, Y, must be dimensioned at 

least equal to NUMPTS in the calling program. 

Controls the number of curves, and/or sets of points 

on one graph: 

This is the only curve, or set of points, to be 

plotted on this graph. 

This is the first of two or more curves, and/or sets 

of points, to be plotted on this graph. 

This is an intermediate curve, or set of points. 

This is the last curve, or set of points, for this graph. 
Controls the type of plot (i.e., curve or point plot): 
This set of points is to be represented by a curve. 

These points are to be plotted with a cross (x). 

These points are to be plotted with a plus (+). 

These points are to be plotted with a square (QO). 

These points are to be plotted with a diamond (0). 


These points are to be plotted with a triangle (A). 


a9 


Dee 


LABEL: 


This is a Hollerith curve or point identifier. 

If a curve is being drawn, LABEL must have 4 char- 

acters (including any blanks), and these will be 

reproduced beside the end of the curve. This 

argument can be set in the calling program by a 

declaration statement, e.g., 

REAL LABEL/4H1234/ 

or REAL LABEL/'1234!'/ 

or REAL LABEL/4H / (The latter must be 
used when no label is 
wanted.) 

If a set of points is being plotted, LABEL is an 8- 

character identifier. The first 4 characters will be 

reproduced beside the first point, and the last four 

characters will be reproduced alongside the last 

point. This argument can be set by a declaration 


Seaeement., E€.e., 


REAL*8 LABEL/8HFRSTLAST/ 

or REAL LABEL(2)/4HFRST , 4HLAST/ 

or REAL*8 LABEL/8H / (The latter 
must be used 
when no label 
is wanted). 

The above arguments, a. through f. (and q.), have 

meaning every time DRAW is called. On the other 

hand, the remaining arguments, g. through p., have no 


meaning except when MODCUR = O or Il. 
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ee 


5 Ga lit ba 


EXSCAL;: 


YSCALE : 


TLTAUP: 


LYRIGH: 


MODXAX : 


An array of 96 EBCDIC characters, the first 48 of 


which will form the first title line, and the last 


48 the second. The array must be dimensioned in 


the calling program, must contain the user's job 
identification, and must have unwanted characters 
set to blank. This may be done with a declaration 
Statement such as: 

REAL*8 ITITLE(12)/'96EBCDIC character'/ 
X-scale in units per inch. EXSCAL will always be 
rounded off to one figure significance. If EXSCAL = 0, 
the X-scale will be computed by DRAW. This is called 
auto-scale. 
Y-scale in units per inch, with properties corresponding 
to those of EXSCAL. 
Distance, in inches, of the X-axis from the bottom of 
the graph (OSIXUP SIHIGH). This argument will be 
ignored unless MODXAX = 2, see below. 
Distance, in inches, of the Y-axis from the left of 
the graph (OSIYRIGHSIWIDE). This will be ignored 
unless MODYAX = 2, see below. 
Determines the mode of the X-axis location: 
The X-axis will be located automatically by DRAW, 
with the origin of Y on the graph. 
The X-axis will be automatically located by DRAW, 
with the origin of Y moved (in one's imagination) an 
integer number of inches above or below the graph, if 
this is appropriate. This option can be used only if 


the Y-scaling is automatic (YSCALE = 0). 
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M 


MODYAX : 


IWIDE: 


IHIGH: 


IGRID: 


LAST: 


The X-axis location will be as specified by IXUP. 
Determines the mode of Y-axis location in the same 

way as MODXAX, above, governs the X-axis location. 
Width of graph in inches (1SIWIDES9). If IWIDE is 
out of this range, a value of 8 will be assumed. 
Height of graph in inches (IS IHIGHS15). If IHIGH 
is out of range, a value of 8 will be assumed. 

If IGRID = 1, a 1" x 1" grid will be superimposed on 
the graph. This is useful only if plain paper is 

used on the CalComp Plotter. 

Indicates to the calling program whether the previous 
plot was completed successfully. The codes are: 

Last plot was completed successfully 

Last plot was not completed successfully. 

Last plot was not completed successfully and no further 
graph output will be attempted until DRAW is next 
entered with MODCUR = 1 or 0 

An attempt was made to enter DRAW With MODCUR # 1 or 0 
while the error lockout was set. This argument must 
always be a variable name and never a number in the 


call statement. 
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oy. 


Notes and Comments: 


ae 


The graph scales and, if MODXAX = 1 and/or MODYAX = 1, the 
amounts of origin offset are always output as part of the 
graph title. 

There are internal checks of the input to DRAW to prevent 
incorrect use. If an input error is detected, an attempt 
will be made, where possible, to complete the plot. If 

an argument is “corrected” in this process, the user will 

be so informed on the standard output. If it is not possible 
to complete the plot, the user will be informed of the 

reason by a message on the standard output. 

If part or all of a curve would fall more than 0.6" laterally 
beyond the ends of the X-axis, or 0.7" vertically beyond 

the ends of the Y-axis, the X and/or Y ordinates will be 
limited so that the curve will typically become a line along 
part or all of the boundary of the graph as‘here defined. 

If one or more points of a point plot would fall outside the 
graph area, the plot of the point, or points will be inhibited. 
The number of such points will be reported to the user on the 
standard output. 

It should be pointed out that the X and Y scaling and axis 
locating processes are entirely independent, so that, for 
example, X might be auto-scaled, while the Y-scale is 
specified. At the same time the X-axis might be located 
automatically, while the Y-axis location is specified. 

It must be remembered that the scales and axis locations of 
a multi-plot graph are set when DRAW is called for the first 


time (with MODCUR = 1). Thus the user must attempt, at the 
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time to achieve scaling and axis location which will be 
appropriate to all the plots he intends to make on the 
one graph. Particularly if the automatic features of 
DRAW are selected, foresight will be demanded of the user 


in this respect. 


Auto-Scale Properties: 

The scale factor is chosen from amongst the values 1, 2, or 5 
units per inch, or some power of 10 times one of those values. 
A curve, or set of points which is plotted with auto-scale will 
normally lie entirely within the graph area is defined in 3.d., 
above. The only exception may occur if an axis is placed, by 
the user, along one edge of the graph (e.g., IXUP = 0, MODXAX = 2). 
In ch a case, points “outside” the axis are not considered in 
the selection of a scale factor (e.g., negative Xx. do not affect 
the choice of scale when IXUP = 0). If automatic axis location 
asevelivas aubu[ceale 15 Selected, Ehewploc, 1f 1t does not fim 
the graph area, will be placed as far as possible towards the 
bottom-left of the graph area consistent with the fact that the 
axes can be set in increments of 1" only. 


lst revision - 31 July 1967. 
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APPENDIX B 
COMPUTER PROGRAM FLOW CHART 


DIMENS ION 
READ FIRST SET OF DATA 


IS E=0 7 = +! Go To END 
Oo 


N 


LOAD PARAMETER 
INITIALIZE 


SUBROUTINE FACTOR 












Select factors 
for secondary ef- 
fect combination 





PRELIMINARY CALCULATIONS 


ASSUME INITIAL VALUES OF Fts 


COMPUTE W AND W 






SUBROUTINE ALPHA 
Compute ots or O's 





COMPUTE NEW VALUES OF F'S 


TEST OLD F'S AGAINST NEW F'S 


IS CONVER-= 
GENCE CRITERIA 
SATISFIED ? 




















IMIPRCVE 
VALUES 


Ole 'S 
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COMPUTE MOMENT ETC, USING 
CONVERGED VALUES OF F'S 





ASSIGN CONVERGED VALUES OF 
F'S BOP EROreR FTE 


IS TABULAR : 
OUTPUT SPECIFIED ? 


9 








. 


PRINT RESULTS 


HAVE ALL COMBI- 
CHANGE No NATIONS OF SECON} 
COMBINATION DARY EFFECTS BEEN 

INVESTIGATED ? 


WES 


HAS ENTIRE RANGE 
OF PARAMETER BEEN )-%2 pre 
INVESTIGATED ? 


YES 


IS GRAPHICAL OUT- 
BUT SPECIFIED ? 
Wes 


CALL DRAW 


NO 
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APPENDIX C 


DATA CARD DESCRIPTION 


C-l1 Data Card Placement 

For one set of data, six data cards are required to give the 
program the necessary input. The first card gives the physical data 
such as leg length, pipe size and material, and other quantities that 
describe the physical situation. The second card supplies the control 
variables to the program. These control variables indicate to the 
program what parameter is to be investigated and its range. 

The last four cards supply the initial values given to the 
FIE's (see page 26). Each of the last four cards contain eight of 
these values. Unless a particular computer run starts near the 
critical part of the parameter range (where convergence is difficult) 
the values on the last four cards are usually zero. If such is the 
case, blank cards may be used; however, when values other than zero 
are supplied as input, the format shown in Appendix C-2 must be used. 

It must be remembered that six data cards must always be used 
for a particular set of input whether or not the FTE's are zero. 

Any number of sets of six data cards can be supplied to the 
program for any one run. Consideration should be given to time 
requirements, however. After the final set of six data cards, one 
blank card follows and then, after that, the last control card (see 
Appendix D). The blank card is used to terminate the run by telling 
the program that all data has been considered. 

The user should be familiar with fixed and floating point variables 


and the common Fortran field specifications. 


3 Daniel D. McCracken, A Guide to Fortran IV Programming (New York: 
John Wiley and Sons, Inc., 1965), pp. 3-5, 13, 84. 
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Appendix C-2 shows a typical input (all six data cards). The 
extreme left hand column indicates the data card number (this order 
must always be followed). For each data card there are two rows 
shown. The top row indicates the variable name whose value is 
shown directly below in the next row. This second row is the 
actual value of the variable named above it and appears just as the 
data card itself does. The top row of the two in each case serves 


only to aid in identification of the variable. 
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C=2 Sample Data Card Set 








le VOMICW) ee ee ere. apa ‘sala tT, alsalegiea el zel afogiestactastacleciesleelz7e ae Ae 2o|99 Solovlcelzel iwpowlscieciselorisetesicgiz2s! efor ezieziezioziezieziczizziizpozie: ‘ets Silgilelgilaisjole! a | 
| * | a fe i | Laie T° }6F ieee (cde Shot Ih lee lias epee 2) cad gt aa Gia a | 7 i PF Sa |e ere fie ata | bet fags os ed 
I | 

| eB 7 ee Sime to PF ot taht Tf 7 i Gt RSis an) RIS, a ia 0d ae See emi ipa vas Seats a ar ae Tic, UL a lem (a Ue SG) ed i 
TT ee Le | 0a i aa Cm ay FR OG a (ee a ie ows lriehorlgale ae Ltr th eer! Ter Pict a 
‘Od Ba 4 roleol Igy CT ior® OW! Tg T aia’e ry Olo' gh TTT 'e ood iqitiya o”' ‘id O 00 

an ; a 

somata | avaaas | esigaa sTalz 2prsig Wo) sera 

tat t 1) vee tao. ery le PF eer ieee <a a leas Da of (ia (cai Sa (as im 


oO TD '' 1 Ofbd GIT O'bd git PE Foo GT TT Flo MITT Epo igettitot6 og ag ittavwo d 
“=] > [aati TT vrgsgly ali "Cli ay' | Wt i ie VWrpsig TTyOrfaly' " ~algosti ‘yorsrt 264!- 
amici. | om bl ier tr td | ee toy Tv tal Te Souk 1 

oo iq TTT ig od WT TTT H*Iloo G7 1) Ig *Ioor arr gtgo grt Toor a. ooo’ 7d 

“Tae ay eae C733" re a er 
Ty ter 1 1 ae ler) te "chomet | | tal oe! alot mis bot 1 je © 

OO dq TTT ojo WdittT oppo grit ojo gtr ©¢ 00" ‘aT ovlou yo boa 
oe a (CMS | era eas OTS, ue 

; T° tile fot tele ig G(R aa cas uO ra 1 i leer roe 

le 1 (Fa SF I fy 2 ° Gale Joy yA ae 9 “le re a7" =] ia i 2 lie if [ all J 


Se etege ee BT | LIA OW We WA NT | | Tevas 





i a se (imimeioe! gel) load 7 aa | a 
| 
o-g ' ez eo GT YLiay Tt oooy AO 1 | OFZ Soy! WMolsoloI 7 11 olsieigigr T 7 > ojro9e | r 
. ‘I4d35gT ae ae - TIT Terr TT@WSHLO 77 Logs TINT @ erry cero TT Stat ArT 









va egos ss Tou leet octee] 721 2 [621591 aol 29! 991 sol oa] colzol ol osleslastesloslssleslecizsl is loslevlavlze lovlsvivelcolzel is lovles les lig loclccl oe lsclzetis 06 lezlarlazlerlsziczlcviraliz az si lalalolelelelalyviotelelilole 





Z WHdly :2 I VHdlv -I 
Y¥38SWNN OML =:2 ano - | 
VWt!tH3s - ee —emme ee 
iNAWw B1vis. NVYLYOS 


O VWHdlV = 























ee ee 


(alae SP 
%ISS3308d YO4 LINENS 


a ee ne, 


14a 3d Td 











31v0 


Wy HDOH 


BLANK CARo —} 








CG 2 

ONE SET vf 5 ———ne 
OF a 2 
SIX 03 a 


DATA 
CAROS 









SUBROUTINE ALPHA ¢ 


MAIN 
PROGRAM 


FORT. SYS IN,A DDA* 


APPENDIX D 


CONTROL CARD ARRANGEMENT 
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E-l Main Program 


APPENDIX E 
PROGRAM LISTING 
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E-3 Subroutine FACTOR 
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APPENDIX F 


SAMPLE CALCULATION 


The following table of values is provided to aid the user in 
checking the computer calculations. Two combinations of secondary 
effects are listed. The first is the "Exact Value" combination 
where all secondary effects are included. The second is the most 
erroneous combination, or the one where all secondary effects are 
suppressed in the calculation, The descriptive data is that listed 
on page 42 of Section IV and also exactly the same as the sample 


data card in Appendix C, with a leg ratio of 0.2. 


Quantity All Effects Allowed All Effects Suppressed 
Assumed Fl 65372.66553 74373.10291 
Assumed F2 4623.212071 5192.148178 

W 0054419338 0.0 

oF 5011360506 5000000000 
oF 3340605064 . 3333933135 
A, 1.001817935 1.000000000 
0Ha 1.002727150 1000000000 
Ot, .4185170188 .4166666865 
HX «2678453147 . 2666666508 
A 0072796731 0.0 

Ag 6695786688 6666666865 
Ay 4320902252 4166666865 
Oro 28594097 26 2666666508 
WBAR 1.923736751 0.0 

esl 2.324235399 .5000000000 
Oe 1.496163752 3333333135 
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Quantity 


Qi XI Q! 


Qtr Qs 
3° ¢ 
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Q OQ 


io 


Fr 
Re 
or 
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Gal 


Jee 
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Fl Converged 
F2 Converged 


Moment 


All Effects Allowed 


3.8/78225194 
5.471216862 
8.613517662 
5 .485230450 
31.81243915 
13.54446691 
8.624985 666 
5 .499880474 
PeOl? Loge 
205.4115482 
- .0001242864 
-0018232014 
0000294904 
- ,0001211686 
-0132454070 
0192106939 
65840.74973 
4633 .905400 


[Oya o0Tso7 1 
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All Effects Suppressed 


1 .000000000 
1.000000000 
-4166666865 
- 2666666508 
Oro 
6666666865 
-41666668 65 
- 2666666508 
- 6000000000 
150.0000000 
- .0000882817 
0013242253 
-0000247189 
- 8828170315 
0.0 
0,0 
74972.17486 
520gn 96999 


12390822 . 149 


APPENDIX G 
SHORTENING DUE TO DIFFERENCE BETWEEN ARC AND CHORD 

The secondary effect produced because of difference between 
arc and chord is described by equation (7) of Section II when 
higher order terms are neglected. The following is a more detailed 
development of that relationship. 

In order to have a physical understanding of the background of 
the derivation, refer to Fig. G-l. Due to the action of cross shear, 
as it is idealized in customary theory, the horizontal length, dx, 
from 0 to A,is stretched as it is inclined to angle YX, so that it goes 
from 0 to B, Then the action of bending further inclines it so that 
it takes the position OC, but there is no change in length during this 
bending process. (Change in length due to axial force is accounted 
for separately.) The increment of difference between arc and chord 
is the distance dd shown in the figure. 


In the derivation, it is first assumed that ¥ is small so that 


tan¥ =¥ approximately. 


ds*= dx + dy? 
dr = dx- ds cos(y) 
zr. = Y or Cly = ¥ dx 
ds*™ = dx? + ¥*adx? 
ds = dxt (1 +¥) 
ass dxvity® 

d= dx- dx iF ¥? Cos(y’) 

dh= dx (1-(¥? cos(y’)) 





ot 


If series representation of cos(\) and Vity? is substituted 


with higher order terms neglected 


dn = dx(1-Crd¥"(i-4b))) = 88 ((y)E Y?) 


(7) N= , l(y) - *) dx 


Or 
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